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SECTION  I 


INTRODUCTION 

The  immediate  purpose  of  calculations  presented  here  is  to  provide 
predictions  of  effects  associated  with  earth  pe.  etration  experiments  con- 
ducted at  the  Watching  Hill  site  in  Alberta,  Canada,  in  mid-July  1974- 
In  particular,  predictions  of  deceleration  history  and  medium  stre  s 
histories  at  specific  locations  were  desired  for  comparison  with  experi- 
mental measurements  of  these  quantities.  Such  comparisons  permit  evalua- 
tion of  current  theoretical  prediction  capabilities  and,  in  addition, 
provide  insights  into  the  primary  mechanisms  and  relevant  phenomena 
governing  projectile  penetration  into  earth  media.  Ultimate  objectives 
of  the  computational  effort  are  (l)  to  develop  methods  for  confidently 
predicting  loads  on  earth  penetrating  projectiles  in  order  to  facilitate 
rational  design  of  penetrators  with  a minimum  reliance  on  empirical 
testing,  (2)  to  develop  an  understanding  of  the  detailed  mechanisms  in- 
volved in  earth  penetration  to  the  extent  that  the  empirical  soil  factor 

1 2 

of  Young's  penobration  equations  ’ maybe  correlated  with  appropriate 
target  medium  properties,  and  (3)  to  establish  the  degree  of  complexity 
in  description  of  target  medium  properties  necessary  for  accurate  calcu- 
lation of  projectile  decelerations,  loads,  and  depths  of  penetration. 

Two  finite  difference  code  calculations  were  made  to  estimate  results 

3 4 

of  the  Watching  Hill  experiments.  The  TOODY  Lagrangian  code  was  employed 
to  predict  penetration  performance  of  both  a rigid  and  a deformable  pro- 
jectile. In  both  cases,  the  target  medium  properties  were  described  by 
using  identical  soil  cap  models.  With  the  projectile  considered  as  a 
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rigid  body,  the  finite  difference  code  is  capable  of  more  efficient 
operations  than  when  the  projectile  is  taken  to  be  a deformable  body. 

Rigid  body  calculations  then  permit  examination  of  penetration  to  deeper 
projectile  depths  in  reasonable  amounts  of  computer  time.  To  verify  the 
validity  of  the  rigid  body  assumption,  a deformable  body  calculation  was 
performed  until  the  nose  of  the  projectile  was  fully  embedded.  From  the 
deformable  body  calculation,  stress  and  strain  histories  and  distributions 
within  the  projectile  are  obtained. 

5 

Projectile  penetration  experiments  at  the  Watching  Hill  site  in- 
cluded one  for  direct  comparison  with  calculations  in  which  the  projectile 

impacted  the  target  medium  at  a velocity  of  152.4  m/s  (500  ft/sec).  The 
<5 

projectile,  fabricated  from  D6AC  steel,  had  a diameter  of  O.165  m (6.5 

inches),  a length  of  1*524  m (60  inches),  and  a mass  of  181.44  kg  (400 

pounds).  Nose  shape  was  9*25  CRH  tangent  ogive  with  a slight  b]unting 

of  the  tip.  Accelerometers  mounted  within  the  projectile  provided 

deceleration  histories  during  penetration.  Integrations  of  deceleration 

histories  result  in  projectile  velocity  and  position  information  including 

depth  of  penetration.  Projectile  deceleration  and  penetration  depth  were 

felt  to  be  insufficient  information  for  critically  evaluating  results  of 

computer  code  calculations,  particularly  the  material  property  models 

7 

employed  in  the  codec.  A much  more  severe  test  of  target  medium  con- 
stitutive models  is  provided  by  measurements  of  the  medium  response.  For 

8 9 

these  reasons,  lithium  niobate  pressure  or  stress  transducers  were 
incorporated  into  the  Watching  Hill  experiments.  Six  stress-history 
transducers,  four  of  which  were  to  measure  mean  stress  and  two  radial 
stress,  were  emplaced  in  the  target  medium  in  two  vertical  holes  near  the 
vicinity  of  the  projectile  penetration  locations.  Positions  of  the  stress 
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transducers  are  Indicated  in  Figure  1-1.  At  the  position  (range/depth 
relative  to  the  projectile  impact  point)  of  0.45  m/l.8  m,  a mean  stress  ^ 
transducer  is  located.  A mean  stress  and  a radial  stress  transducer  are 
located  aJ-  0.45  m/3. 6 m,  and  at  0.45  m/5.4  m,  two  mean  stress  gauges  are 
placed.  A single  radial  stress  gauge  is  located  at  0.9  m/3.6  m.  Gauges 
shown,  in  Figure  1-1  at  a range  of  1.5  m were  used  for  another  projectile 
penetration  experiment  but  are  physically  the  same  set  of  gauges  shown 
at  a range  of  0.45  m le  experiment  of  primary  interest.  Measurements 

of  target  response  a. - 

before  been  male;  the  successful  recording  of  medium  stress  histories" 
during  the  Watching  Ilill  penetration  experiments  has  provided  valuable 
information  for  more  direct  verification  of  computer  code  calculations. 

Output  desired  from  computer  code  calculations  has  been  specified 
11 

in  detail  and  is  briefly  summarized  in  Figures  1-1  and  1-2  for  the 
rigid  body  calculations.  Similarly,  for  the  deformable  projectile  cal- 
culations, desired  output  is  indicated  in  Figures  1-3  and  1-4. 


.jectile  penetration  apparently  have  never 
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are  required. 
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SECTION  II 

MATERIAL  MODELS  FOR  TARGET  MEDIUM 

The  Watching  Hill  Blast  Range  of  the  Canadian  Defense  Research  Estab- 

5 

lishment,  on  which  site  the  earth  penetrator  experiments  were  conducted, 

is  located  about  30  miles  north  of  Medicine  Hat  near  Suf field,  Alberta. 

The  target  medium  at  Watching  Kill  consists  of  a thick  succession  of 

glacial  lake  deposits,  outwash  materials,  and  glacial  tills;  the  near 

surface  deposit  is  composed  of  thin  interbedded  layers  of  lacustrine  silt, 

12 

sand,  and  clay  sediments.  The  medium  is  porous  and  contains  varying 

14 

amounts  of  water,  increasing  with  depth.  For  purposes  of  modeling  the 
target  medium  in  computer  calculations  of  projectile  penetration,  the 
medium  has  been  represented  by  four  uniform,  homogeneous  isotropic  layers. 
The  three  uppermost  layers  were  each  taken  to  be  2.4  m (7.87  feet)  thick. 
The  fourth  layer,  beginning  at  7-2  m (23.6  feet),  is  below  the  water  table. 
Material  properties  derived  from  uniaxial  strain  and  triaxial  stress 
experiments  have  been  specified  '3s  for  each  of  the  four  layers. 


■ft 

Of  course,  many  more,  and  much  thinner  layers  actually  exist  in  the  region 
of  interest.  It  was  felt  that  the  short  time  available  for  material 
modeling  made  it  necessary  to  average  properties  into  the  sturcture  chosen, 
using  overall  similarities  of  adjacent  materials  as  guides.  It  w as  also 
necessary  in  some  cases  to  estimate  the  effects  of  loading  rates  on 
material  response.  Therefore,  the  term  "data"  appearing  in  this  section 
should  be  taken  to  mean  what  is  commonly  referred  to  as  "recommended 
properties ." 


For  the  calculations  of  projectile  penetration  with  both  a rigid 

and  deformable  projectile,  a soil  cap  model  is  used  to  describe  each  of 
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the  four  layers  of  the  Watching  Hill  medium.  The  cap  models  employ  a 
failure  envelope  or  fixed  plastic  yield  surface  represented  by  a function 
of  the  form  f^(T,p)  = 0 in  which  p is  the  mean  stress  or  pressure  and  t 
is  the  octahedral  shear  stress  (see  Figure  II-l) . Octahedral  shear  stress 

i 

is  defined  by  t = (o',  .o'.  ./3)2,  where  the  deviator  stress,  o'.  . = a.  . - 
p6^.  The  "cap"  on  this  surface  is  given  by  f^p,^^)  = 0,  where  is 
the  plastic  volumetric  strain.  Position  of  the  cap  on  the  failure  envel- 
is  determined  by  e^,  and  the  cap  expands  or  contracts  as  increases  or 
decreases.  Plastic  strain  increments  are  determined  from  f^  and  f by 
means  of  the  associated  flow  rule.  Figure  II-2  illustrates  an  example 


Figure  II-l.  Cap  Model  Illustrating  Failure  Envelope, 

Cap,  a Stress  Path,  and  Cap  Contraction 

of  the  hysteresis  produced  by  a cap  model  description  in  a series  of  load- 

unload  cycles  of  a uniaxial  strain  test  and  the  corresponding  path  in  t - 

12  ... 
space.  Information  supplied  for  the  Watching  Hill  target  medium  consis- 


ted  of  a single  uniaxial  strain,  load-unload  path  and  a failure  envelope 
considered  to  be  representative  of  each  of  four  layers  of  the  medium. 

From  the  data  available  for  each  of  the  idealized  four  layers,  parameters 
were  determined  for  the  functions  f ^ and  fg  . For  Layers  1,  2,  and  3?  "the 
failure  envelope  is  given  by 

t = A - C exp(-Bp) 

and  for  Layer  4 by 


t = A[l  - 3p/b)2]  + C,p  < B/3 
t = A + C,  P 2 B/3 

where  p is  pressure.  The  yield  surface  ellipse  for  Layers  1,  2 and  3 is 
described  by 

T = [(X  - l)2  - (Vbp  - l)2]Vr  > 

where  X represents  the  intersection  of  the  cap  with  the  p axis  and  is 
given  by 

X = -L-Vn(l  - 6p/w) l/k 

L represents  the  value  of  p at  which  the  cap  and  failure  envelope  intersect, 
and  at  which  the  cap  has  a horizontal  tangent  and  is  given  by 


(L  for  1^0 

L(0  = 

p (0  for  l < 0 

where  l is  the  solution  of  the  transcendental  equation 


l + R{A  - C exp(-Bl/^~3)}  - X(sp) 
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Figure  II -2.  Loading  and  Unloading  Cycles  for  Cap  Model 
Material  in  Stress-Strain  and  Shear  Stress- 
Pressure  Space 


For  Layer  4,  the  yield  surface  is  taken  as 


f = - P/V~2  • 


In  the  calculation,  the  cap  determines  the  plastic  strain  increment 
whenever  p S:  L/^~.  This,  coupled  with  the  fact  that  the  normal  to  the 
cap  never  has  a negative  spherical  component,  precludes  bulking,  or  di la- 
tency, when  the  loading  state  is  on  the  cap. 

The  initial  value  of  the  quantity  X for  each  material  was  taken  to 
be  one  percent  of  the  value  necessary  for  the  cap-envelope  intersection 
point  to  be  at  p = 0.  This  results  in  a very  small  initial  elastic  range. 
As  plastic  deiormation  occurs,  the  growing  cap  remains  centered  on  the 
origin,  until  X exceeds  100  times  its  initial  value. 
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Elastic  bulk  moduli  are  given  by 


K = Bq{1  - Bl  exp(-B2p)l 

for  Layers  1,,  2,  and  3?  and  by 

K = Min{ Bq  exp(B2p/B1);  BQ  exp(B2/VT)} 

for  Layer  4.  The  elastic  shear  moduli  are 

3(1  - 2„)ll  + S1  exP<-s2P)l  „ 

2(1  + v)  Cl  + "s^  * 9 

in  which  v is  Poisson's  ratio.  Constants  for  the  various  expressions  are 

listed  in  Table  II-l.  Figures  II-3  through  11-6  compare  data  obtained 

from  the  cap  models  (which  are  used  in  code  calculations  of  projectile 

penetration)  witn  laboratory  data  from  uniaxial  strain  and  triaxial  stress 

experiments.  Figures  I 1-7  through  11-10  show  data  for  uniaxial  strain 

paths  in  t - p space  and  corresponding  cap  model  results.  The  cap  model 

for  Layer  4 fits  the  data  well;  however,  for  Layers  1,  2,  and  3?  less 

satisfactory  agreement  between  model  predictions  and  experiment  is  obtained 

with  the  model  response  being  less  stiff  in  shear  than  indicated  by  da  "a. 

This  suggests  the  need,  in  matching  static  laboratory  data,  ior  a more 

elaborate  model  than  the  12-parameter  model  used.  One  such  model  has 

been  suggested1^  in  which  the  strain  hardening  cap  is  described  by  an 

ellipse  of  variable  eccentricity.  Cap  model  formulations  with  29  para- 

1.3 

meters  available  for  fitting  data  have  been  described  ’ these  formu- 
lations undoubtedly  can  better  fit  the  available  laboratory  data  and  could 
be  employed  in  future  calculations  to  assess  sensitivity  of  calculated 
projectile  penetration  results  to  cap  model  expressions  and  parr-uete^ . 
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TABLE  II -1 

Cap  Model  Constants  if  r tte  Four  Layers  of  the  Watching  Hill  Medium 


Constant 

Units 

Layer  i 

Layer  2 

Layer  3 

Layer  4 

Po 

(Mg/m3) 

1.490 

1.426 

1.859 

1.971 

Bo 

(GPa) 

O.6237 

O.8889 

1.214 

1.103 

V 

0.31 

0,l4 

0.257 

0.48 

B1 

0 

0.7622 

0.8622 

— 

B1 

(MPa) 

— 

— 

— 

2.382 

B2 

(MPa)-1 

0 

0.2087 

0.3965 

— 

B2 

— 

— 

— 

2.139 

S1 

(MPa)"1 

0 

0.35 

-0.8 

0.51 

S2 

0 

0.3859 

0.3464 

0.6603 

A 

(MPa) 

1.966 

7.6o4 

0.7084 

8.123  x 10' 

B 

(GPa)"1 

346.4 

91.42 

468.5 

— 

B 

(MPa) 

— 

— 

— 

3.102 

C 

(MPa) 

1.902 

7.546 

0.6368 

8.123  X 10' 

D 

(GPa)-1 

38.61 

131.2 

150.0 

800.0 

w 

0.37 

0.16 

0.12 

-4 

1.0  X 10 

R 

2.252 

2.i79 

2.785 
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Figure  II- 6.  Cap  Model  Fit  to  Uniaxial  Strain  and  Triaxial  Failure 
Envelope  Data  for  Layer  4 at  Watching  Ilill 
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Figure  II-7.  Uniaxial  Strain  Path  in  T-p  Space 
for  Layer  1 
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Figure  II- 8.  Uniaxial  Strain  Path  in  j-p  Space 
for  Layer  2 
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It  should  be  pointed  out  that  the  cap  model,  as  implemented  in  TOODY 

for  the  calculations  reported  here,  differs  in  some  details  from  the  cap 

15 

model  described  by  Nelson  et  al.  In  particular,  the  relations  for  the 
variation  of  bulk  and  shear  moduli  with  p and  t do  not  satisfy  the  com- 
patibility requirements  which  arise  from  the  assumption  of  the  existence 

of  an  elastic  strain  energy  function.  Also,  the  particular  function  chosen 
to  model  the  cap  for  Layer  t exhibits  a "corner"  on  the  pressure  axis, 

reminiscent  of  those  appearing  in  the  Tresca  yield  condition,  for  example. 
The  effect  of  these  variations  on  the  particular  calculations  reported 
here  is  difficult  to  assess,  particularly  as  no  purely  spherical  stress 
states  were  encountered  in  the  calculation. 


SECTION  III 

RIGID  PROJECTILE  CALCULATIONS 
3 

The  TOODY  Lagrangian  wave  propagation  code  was  used  to  predict  tar- 
get and  penetrator  motions  and  stresses  during  the  experiments  planned 
for  the  Watching  Hill  test  site.  The  code  is  termed  Lagrangian  because 
the  nodes  of  the  calculational  grid  move  in  response  to  the  local  stress 
gradients.  The  basic  process  the  code  performs  is  to  solve  difference 
analogs  to  the  partial  differential  equations  describing  the  balance  of 
momentum,  mass,  and  energy,  and  the  constitutive  relations.  This  procedure 
is  executed  at  each  node  of  the  grid,  the  time  is  then  advanced,  and  the 
whole  process  repeated  until  the  desired  final  time  of  the  calculation 
is  reached. 

Since  the  time-step  with  which  the  calculation  advances  is  control- 
led by  the  minimum  distance  on  a mesh  over  the  entire  grid,  large  distor- 
tions can  result  in  an  inefficient  calculation.  Large  distortions  can 
also  degrade  the  accuracy  of  the  numerical  methods  used  in  the  code, 

4 

which  then  yield  unrealistic  results.  For  these  reasons,  the  TOORE/i 
code  was  developed  to  permit  carrying  Lagrangian  code  calculations  past 
the  point  where  they  would  otherwise  have  to  be  abandoned.  TOOREZ  allows 
the  user  to  redistribute  the  nodes  of  the  calculational  grid  (i.e.  rezone 
the  problem),  then  repartitions  the  mass,  momentum,  and  energy  In  such  a 
way  that  all  these  quantities  are  conserved.  Using  the  output  from  the 
rezoning  code,  the  calculation  can  then  be  resumed. 

The  positions  of  nodes  in  the  calculational  grid  which  originally 
lie  on  the  axis  of  rotational  symmetry  are  prevented  from  occupying  the 


2b 


r>  «? »&r  'vnwwptf^rjf&m gp^yj^;arg»qs^^ 


space  occupied  by  the  penetrator.  The  motion  of  each  of  these  nodes  is 
first  calculated  in  the  standaru  way  from  the  stress  divergence,  including 
the  effects  of  soil-penetrator  friction,  if  any.  In  this  calculation,  a 
condition  of  perfect  slip  was  assumed  for  the  soil-penetrator  interface. 

If  the  standard  motion  calculation  for  the  node  yields  a position  on  the 
surface  of,  or  outside,  the  projectile,  no  adjustments  are  necessary. 

When  a node  is  calculated  to  lie  inside  the  penetrator,  adjustments  to 
the  velocity  and  position  are  made  so  that  the  normal  velocity  relative  to 
the  penetrator  is  zero,  the  tangential  velocity  has  that  value  already 
calculated,  and  the  position  is  on  the  penetrator  surface.  In  resolving 
the  velocity  into  normal  and  tangential  components  on  the  ogival  nose,  the 
position  of  the  node  at  the  previous  time  step  is  used  to  evaluate  the 
local  normal  direction.  The  areas  and  volumes  for  material  and  momentum 
zones  adjacent  to  the  penetrator  are  adjusted,  where  necessary,  to  include 
the  effect  of  the  penetrator  surface. 

At  the  beginning  of  each  time  step  of  the  calculation,  surface  trac- 
tions for  all  soil  zones  in  contact  with  the  penetrator  are  integrated 
to  give  the  total  axial  force  acting  to  decelerate  the  penetrator.  For 
a zone  in  contact  with  the  penetrator  surface  the  stress  is  assumed  con- 
stant, and  the  axial  contributions  of  the  normal  stress  and  frictional 
stress,  if  any,  are  expressed  as  functions  of  position  on  the  penetrator 
surface.  These  are  integrated  over  the  range  covered  by  the  zone,  so  that, 
on  the  nose,  changes  in  the  local  normal  direction  and  surface  area  are 
accounted  for. 

For  this  calculation,  it  was  desired  to  model  the  target  to  a large 
depth.  Therefore,  the  initial  zoning  of  the  target  region  utilized  expon- 
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tially  increasing  zone  sizes  to  allow  good  resolution  near  the  surface 
and  maintain  a reasonable  number  of  zones  in  the  problem.  Table  ITI-1 
and  Figure  II1-1  describe  the  initial  zoning.  As  penetration  depth 
increases  to  the  legions  where  the  zone  sizes  become  too  large  for  good 
resolution,  the  rezoning  code  is  used  to  adjust  the  zoning,  providing 
smaller  zones  in  front  of  and  on  the  nose  of  the  penetrator.  It  was 
determined  early  in  the  calculation  that  no  significant  changes  occur  in 
the  soil  once  the  nose  of  the  penetrator  has  passed  by,  so  nothing  is 
lost  by  decreasing  the  resolution  behind  the  penetrator. 


TABLE  III-l 

Initial  Zoning  for  Rigid  Projectile  Calculation 

FOR  ALL  REGIONS I z = z.  + ^ (expf(j_  - j)  Ln(l  + \ ))  - l) 

Kz  z 

x - x + ^ (exp( (i  - i ) <,n(l  + \ ))  - l) 

Kx  O .C 

(rn)  X (m)  (m)  Ax(ni)  \ \ j i 

0 q ti  Jt  V 


o.  0.  ,03  .03  0,02  0,02  91  1 

2 

(IS  i < 25; 

l < j < 67)  0.9127  0,  0.0483  .03  .05  .02  67  1 


(26  < i < 46; 

1 < j < 68)  0.  0.91270.04830.0483  0.05  0.05  68  26 


REGION: 


1 

(1<  i<  25; 
67  < j < 91) 


4WM±-**.-  


The  calculation  was  carried  out  to  a time  of  approximately  52.5  ms, 
at  which  time  the  penetrator  had  reached  a depth  of  7. 3 m,  and  retained 
a velocity  of  129  m/sec.  The  peak  deceleration  encountered  by  the  penetra- 
tor was  approximately  55  g's,  and  at  the  end  of  the  calculation  the  deceler- 
ation appeared  to  be  decreasing,  as  the  nose  begins  to  enter  layer  4,  to  a 
value  less  than  20  g's. 

Taking  a constant  deceleration  of  20  g's  in  the  assumed  uniform  materi- 
al of  layer  lj  below  7*2  m leads  to  a value  of  about  0,7  sec  for  total  time 
to  penetration  depth  and  a value  of  50  m for  final  depth  of  penetration. 
These  are  the  results  expected,  based  on  calculated  projectile  motion 
through  the  first  three  layers,  had  the  calculation  been  carried  out  until 
the  penetrator  comes  to  rest.  The  duration  of  projectile  penetration  and 
final  depth  thus  predicted  are  much  larger  than  observed  in  experiment, 
suggesting  a shortcoming  in  the  calculation  due  possibly  to  friction 
effects  or  to  use  of  a material  response  model  for  Layer  4 which  is  not 
an  accurate  representation  of  the  actual  medium  encountered  by  the  penetra- 
tor (see  Section  V for  further  discussion). 

The  calculations  indicated  that  the  value  of  deceleration  achieved 
was  somewhat  sensitive  to  zone  s^ze,  although  the  precise  degree  of  this 
sensitivity  was  not  known.  For  this  reason,  an  analysis  of  the  effect  of 
zone  size  was  undertaken.  Three  calculations  were  run  with  a projectile 
penetrating  0.6  m of  Layer  1,  with  initially  square  zones  of  0.015,  0.03, 
and  0.06  m.  Initial  conditions  were  identical  to  those  of  the  main  problem. 
Figures  III-2  and  III- 3 show  deceleration-depth  curves  and  normal  stress 
distribution  along  the  nose  at  a penetrator  depth  of  0.4  m,  respectively. 

It  may  be  noted  from  tnese  graphs  that  while  alterations  in  initial  zone 
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size  have  a marked  effect  on  the  normal  stress  distribution,  the  decelera- 
tion is  much  less  affected  for  the  range  of  zone  sizes  considered.  This 
is  evidently  rue  to  the  fact  that  where  disparities  in  the  normal  stress 
distributions  are  largest,  the  area  over  which  the  normal  stress  acts  is 
relatively  small.  That  is  to  say  that  for  the  smaller  mesh  calculations 
the  normal  stress  distribution  along  the  projectile  nose  has  its  maximum, 
shifted  toward  the  nose  tip  (see  Figure  III-3) . Figure  III -4  shows  pres- 
sure vs  radial  position  at  the  depth  of  peak  normal  stress,  at  a penetra- 
tor  depth  of  0.4  m.  As  expected,  the  finer  zoning  allows  resolution  of 
steeper  gradients.  Some  sort  of  extrapolation  to  obtain  the  stress  at 
the  soil-projectile  interface,  rather  than  one  half-zone  away,  would 
have  yielded  somewhat  higher  decelerations,  but  pr  bably  not  enough 
higher,  taken  by  itself  to  cause  a large  change  in  duration  of  the  cal- 
culation or  to  influence  significantly  the  values  of  projectile  decelera- 
tion calculated.  A combination  of  considerably  finer  zoning  than  was 
used  an  the  main  calculation,  and  extrapolation  to  find  stress  at  the 
soil- projectile  interface  would  provide  a slightly  more  realistic  cal- 
culation, but  would  also  be  quite  expensive  in  terms  of  computer  time. 
Increases  in  calculated  normal  stress  and  projectile  decelerations  due 
to  finer  zoning  and  extrapolation  are  estimated  to  amount  to  at  most  2Oaj0 
using  the  techniques  employed  in  the  rigid  body  calculation. 

It  was  hoped  that  the  rezones  of  the  problem  could  be  accomplished 
with  no  significant  effects  on  the  histories  being  accumulated  at  the 
selected  points  in  the  soil.  This  goal  is  really  only  a reasonable  one 
for  those  quantities  and  their  integrals  which  are  conserved  during  the 
rezoning  process.  The  higher  degree  of  smoothness  of  the  velocity  his- 
tories as  compared  to  the  stress  measure  histories  reflects  this  fact. 
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Rezones  of  the  problem  were  performed  at 

TIME  (ms)  DEPTH  (m) 


1.38 

.21 

3.05 

.46 

8.0 

1.2 

11.9 

1.78 

l4.6 

2.18 

23.4 

3.45 

33.8 

4.88 

4l.6 

5.91 

Due  to  tape  handling  problems,  history  information  for  the  period 
from  11.9  His  to  14.6  ms  was  lost. 

Appendix  A presents  results  of  the  calculation  in  graphical  form. 
Figures  A-l  through  A-5  show  deceleration,  velocity  and  depth  vs.  time, 
and  deceleration  velocity  vs.  depth.  Figures  A-6  through  A-9  show  pro- 
jectile kinetic  energy,  soil  kinetic  energy  due  to  radial  velocity,  soil 
kinetic  energy  due  to  vertical  velocity,  and  total  kinetic  energy  of  the 
soil,  all  per  radian,  and  all  vs.  time.  Table  A-l  and  Figure  A-10  pro- 
vide keys  to  figures  A-ll  through  A-78,  which  show  histories  of  various 
quantities  at  selected  points  in  the  soil.  (Where  they  appear,  radial 
and  axial  stress  components  are  plotted  positive  in  compression) . Fig- 
ures A-79  through  A-l4l  show  the  calculations!  grid  near  the  projectile, 
normal  and  shear  stress  distributions  along  the  projectile,  and  contours 
of  selected  quantities  in  the  soil.  Due  to  the  rezoning  method  used,  only 
the  grid  lines  with  depths  near  that  of  the  nose  region  are  representa- 
tive of  target  deformation.  Also,  for  projectile  depths  greater  than 
0.4  m,  the  spatial  scaling  is  distorted,  so  that  radial  distances  appear 
twice  as  large  as  their  actual  value.  The  grid  plots  contain  "x"  in  the 


center  of  a zone  calculated  as  failing  in  tension,  and  " c"  where  the 
stress  state  is  on  the  failure  envelope.  When  a mesh  is  calculated  to 
fail  in  hydrostatic  tension,  its  pressure  is  set  to  zero,  and  its  stress 
deviator  is  limited  by  the  smaller  of  the  deviator  stress  on  the  cap  and 
the  failure  envelope  at  :zero  pressure.  Therefore,  a zone  may,  in  some 
cases,  contain  both  x and  c.  It  should  als~  bo  noted  that  the  normal  and 
shear  stresses  plotted  are  actually  those  at  zone  centers;  i.e.,  one-half 
zone  away  from  the  projectile.  For  these  plots,  a positive  normal  stress 
and  a negative  shear  stress  would  retard  the  projectile. 

Table  A-2  provides  a key  to  Figures  A-79  through  A-14ls  It  should 
be  noted  that  the  sign  convention  for  stress  components  is  opposite  that 
used  for  the  history  plots;  i.e.  consistent  with  a positive  normal  stress 
being  one  of  extension. 


SECTION  IV 
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DEFORMABLE  PROJECTILE  CALCULATION 

To  test  the  assumption  of  a rigid,  body  and  to  calculate  the  response 
of  the  vehicle  to  the  imposed  loads,  a second  TOODY  calculation  was  per- 
formed. The  vehicle  was  zoned  as  shown  in  Figure  IV-1,  and  the  ground 
was  defined  by  30-mm  squares  to  a radius  of  0.6  m and  a depth  of  0.9  m. 
For  this  calculation  the  U50  conical  tip  was  included,  as  on  the  actual 
vehicle,  rather  than  the  idealized  perfect  ogive  used  in  the  rigid  body 
calculation.  Figure  IV-1  also  shows  points  in  the  projectile  at  which 
the  progress  of  the  calculation  was  recorded  and  plotted.  Plots  were 
also  made  for  points  in  the  soil  corresponding  to  depths  of  0.1  rn,  0.2  m, 
and  0.4  m,  with  distances  from  the  central  axis  of  75  mm,  O.lU  m,  0.28  m, 
0.44  m,  and  O.58  m. 

The  spacer  was  treated  as  an  aluminum  shell  containing  a vacuum. 

The  electronics  package  was  modeled  with  a hydrodynamic  "foam"  equation 
of  state.  The  rest  of  the  vehicle  was  treated  as  steel.  The  ground  was 
represented  by  the  cap  model  for  Layer  1.  The  input  parameters  for  the 


vehicle  materials  were 
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The  last  four  parameters  pertain  to  the  p-a  foam  model  ; cg  is  the 
ambient  value  of  the  sound  speed  in  the  virgin  foam  at  the  initial  disten- 
tion an>  p is  the  elastic  limit,  and  p is  the  compaction  pressure  for 
which  all  voids  are  assumed  to  be  closed.  Therefore,  c^  is  the  bulk 
sound  speed  of  the  fully  compacted  solid. 

Friction  at  the  vehicle-soil  interface  was  neglected,  and  separation 
of  soil  from  projectile  surface  was  not  permitted.  The  calculation  was 
run  to  a time  of  2,4-0  msec,  at  which  time  the  nose  tip  had  reached  n depth 
of  0.398  m.  This  calculation  consumed  2.94  hours  of  CDC-6600  central 
processor  time.  The  problem  was  not  rezoned. 

Figures  B-l  through  B-7  in  Appendix  B describe  the  configuration, 
velocities  and  stresses  0.40  msec  after  impact,  with  the  nose  tip  0.094 
m below  the  ground  surface.  The  situation  at  1.10  msec  is  shown  in 
Figures  B-8  through  B-l4,  corresponding  to  a penetration  depth  of  0.20  in. 
Figures  B-15  through  B-21  give  the  final  configuration.  Velocity  and 
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stress  contours  are  not  shown  for  the  projectile  material  because  the  high- 
frequency  ringing  of  the  vehicle  renders  contours  there  incomprehensible. 

A few  contours  are  shown  as  passing  through  the  projectile  region;  these 
are  intended  to  clarify  the  relationship  of  the  contours  and  are  unrelated 
to  actual  values  in  the  projectile  material. 

Figures  B-22  through  B-30  display  motion  history  for  three  points 
in  the  projectile.  The  motion  of  points  away  from  the  tip  is  dominated 
by  ringing  that  masks  the  rigid  body  deceleration  of  the  vehicle.  Point 
5,  40.5  mm  back  from  the  tip  and  6.35  mm  from  the  axis,  experienced  an 
initial  acceleration  associated  wit  a deformation  of  the  nose  on  impact, 


followed  by  a deceleration  as  the  nose  absorbed  the  initial  ground  resis- 
tance (Figure  B-28) . 

Stress  histories  at  selected  points  within  the  vehicle  are  shown  in 
Figures  B-31  through  B-36.  In  these  and  subsequent  plots,  "J2"  is  the 
label  for  the  quantity  yR  .o'.  which  was  the  requested  output  quan- 
tity?  where  o'  is  tne  deviator  stress  tensor*  This  quantity  is  propor- 
tional  to  octahedral  shear  stress,  \Jo'.  ,ct'  / 3.  At  a point  in  the  projec- 
tile  6l  mm  from  the  tip  and  12  mm  from  the  axis,  the  peak  mean  stress 
exceeded  0.14  GPa  and  ^4  a .0  [ . exceeded  0.28  GPa  (Figures  B-35  and  B-36). 
At  the  tip  of  the  vehicle  (not  plotted),  mean  stresses  in  excess  of  0.5 
GPa  and  axial  stresses  in  excess  of  1.3  GPa  were  noted  both  in  compres- 


sion and  tension,  and  plastic  yielding  occurred. 

Figures  B-37  through  3-39  contain  normal  stress  histories  at  three 
points  along  the  nose  of  the  prijectile  (see  Figure  IV-l) ; the  shear 
stresses  were  assumed  to  be  zero.  Because  stresses  were  not  calculated 
at  the  interface,  the  values  given  were  calculated  within  the  vehicle, 
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one-half  mesh  width  from  the  surface.  The  largest  stress  in  these  fig- 
ures is  43  MPa,  which  occurred  at  the  point  closest  to  the  nose  tip 
(Figure  B-37) . 

The  total  kinetic  energy  calculated  for  the  soil  followed  the  his- 
tory shown  in  Figure  B-40.  The  saw-tooth  pattern  is  caused  by  numerical 
noise,  a peak  occurring  each  time  the  nose  tip  encounters  a new  grid 
point  in  the  soil.  The  noise  is  accentuated  because  the  total  kinetic 
energy  in  the  soil  is  small,  because  it  is  concentrated  in  a few  zones 
near  the  axis  of  symmetry,  and  because  the  approximation  used  to  sum 
the  kinetic  energy  has  undesirable  characteristics  when  applied  near  the 
axis  of  symmetry.  This  effect  was  not  seen  in  the  rigid  body  calculation 
(Fig.  A— 9)  because  the  time  interval  chosen  for  recording  the  energy 
happened  to  be  almost  exactly  one  period  for  this  oscillation,  and 
because  the  omission  of  the  conical  tip  allowed  a smoother  solution  in 
the  rigid  case. 

The  history  of  motion  and  stresses  at  the  selected  points  in  the 
soil  are  shown  in  Figures  B-4l  through  B-112.  Histories  were  also 
recorded  at  three  points  with  range  0.6  m from  the  impact  line,  but  the 
mean  stress  never  exceeded  10  kPa,  the  level  required  to  activate  cal- 
culation of  motion.  Before  comparing  these  results  to  the  rigid  pro- 
jectile calculation,  the  times  shown  must  be  increased  by  0.4  msec  to 
account  for  the  nose  tip  being  26  mm  shorter  and  for  the  tip  being  buried 
33  mm  at  the  beginning  of  this  calculation.  With  this  adjustment,  the 
soil  histories  are  in  excellent  agreement.  The  peak  mean  stress  observed 
at  the  closest  point  (range  .105  m,  depth  75  mm)  was  2 MPa,  and  the 
largest  radial  velocity  was  l6  m/s  (Figures  B-4l  through  B-46) . 
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The  total  momentum  in  the  vehicle  was  summed  at  0.2-ms  intervals  and 
then  differenced  to  obtain  the  deceleration  history  shown  in  Figure  IV-2. 

The  deceleration  is  up  to  thirty  six  percent  larger  than  the  results  for  the 
rigid  body  calculation  (Figure  III-2).  This  difference  is  related, 
in  part,  to  the  use  of  an  ideal  ogive  in  the  rig  d body  calculation, 
compered  with  the  45°  tip  appearing  on  the  actual  vehicle  and  included 
in  the  deformable  body  calculation.  The  soil  velocity  reaches  180  m/s 
on  the  45°  nose  tip,  whereas  in  the  rigid  body  calculation  the  soil 
velocity  does  not  exceed  60  m/s.  In  the  case  of  the  perfect  ogive,  the 
soil  must  reach  a velocity  of  at  least  50  m/s  to  get  out  of  the  path  of 
the  oncoming  vehicle.  The  higher  velocity  produces  larger  stresses  in 
the  soil;  with  the  d.aformable  body  stress  components  larger  than  12  MPa 
were  calculated,  whereas  no  stress  in  excess  of  3.5  MPa  was  noted  for 
the  rigid  body  calculation. 

Unfortunately,  there  are  no  hard  experimental  data  with  which  to 
compare  the  calculated  transient  vehicle  response.  However,  the 
experience  with  this  deformable  body  calculation  demonstrates  that  cal- 
culation of  the  early-time  vehicle  loads  and  response  is  a reasonable 
and  straight-forward  application  of  existing  numerical  methods,  in 
sharp  contrast  to  the  difficulties  in  calculating  penetration  depth  or 
to  the  impossibility  of  obtaining  this  information  from  force  law  or 
cavity  expansion  theories. 
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SECTION  V 


COMPARISON  OF  PREDICTIONS  WITH  EXPERIMENTS 

5 

More  than  a dozen  projectile  penetration  experiments  were  conducted 
at  the  Watching  Hill  site  with  projectiles  of  various  diameters,  weights, 
lengths  and  nose  shapes  and  with  projectile  impact  velocities  ranging 
from  93  to  202  m/s  (305  to  660  ft/sec).  Test  number  6 had  an  impact 
velocity  of  150  ra/s  (493  ft/sec),  a O.165I  m(6.5  in)  diameter  projectile 
of  mass  l8l.44  leg  (400  lb)  with  9*25  CRH  tangent  ogive  nose.  This  test 
closely  matches  conditions  employed  in  calculations  and  provided  pro- 
jectile motion  data  for  comparison  with  predictions.  Test  number  4 was 
performed  with  a 0.1524  m (6  in)  diameter  projectile  of  mass  127.12  kg 
(280  lb),  6.25  CRH  tangent  ogive  nose  and  impacted  at  a velocity  of  158.5 
m/s  (520  ft/sec).  From  this  test,  medium  stress  histories  were  obtained 
with  which  calculated  results  can  be  approximately  compared. 


Projectile  Motion 

5 

Projectile  motion  measurements  were  obtained  from  accelerometers 

mounted  within  the  projectile  and  consist  of  deceleration  versus  time 

during  penetration,  together  with  projectile  velocity  and  position  his- 

5 

tories  derived  from  integrations  of  deceleration  data.  Medium  response 

10  9 
measurements  were  made  with  six  lithium  mobate  stress  transducers 

at  four  locations  within  the  top  three  designated  layers  of  the  target 

medium.  From  projectile  deceleration  data,  information  is  obtained 

regarding  the  nature  of  the  resistive  forces  offered  by  the  soil  medium 

in  retarding  the  projectile's  motion.  Also,  together  with  calculations, 

inferences  regarding  frictional  forces  acting  on  the  projectile  may  be 


made . 


Further,  from  the  deceleration  history  indications  are  obtained  of  the 
gross  layering  structure  within  the  target  medium.  From  dynamic  stress 
data  amplitudes,  durations  and  arrival  times  are  obtained  of  the  pressure 
disturbance  produced  in  the  medium  by  projectile  penetration.  The  com- 
parison of  this  information  with  calculated  results  provides  one  measure 
of  how  accurately  the  medium  response  is  modeled  in  calculations. 

Figures  V-l  and  V-2  illustrate  projectile  decelerations,  calculated 
and  measured  (Test  No.  6,  real  time,  200  Hz  LPF) . The  rigid  body  cal- 
culation with  zero  friction  predicts  a deceleration  of  about  55  g's  which 
remains  roughly  constant  during  projectile  penetration  through  layers  1 
and  2 (each  2.4  m thick).  After  projectile  entry  into  layer  3 (also  2.4 
m thick)  calculated  deceleration  falls  to  abow  ^5  g's  and  in  layer  4 
deceleration  again  drops  to  20  g's  or  less.  This  stair-step  result  is 
related  to  the  material  properties  assigned  to  each  of  the  four  idealized 
layers  (Section  IX) . Layers  1 and  2 have  shear  strengths  wnich  differ 
at  most  by  a factor  of  2 at  pressures  below  10  MPa  (see  Figures  II-3  and 
II-4) . Layers  3 and  4 both  have  shear  strengths  much  less  than  layers 
1 and  2 and  the  strength  of  layer  4 material  is  considerably  less  than 
that  of  layer  3 over  the  pressure  range  of  interest  (see  Figures  II-5 
and  II-6) . Assuming  that  the  soil  medium's  resistance  to  shear  deforma- 
tion is  the  primary  mechanism  influencing  projectile  deceleration  level 
then  the  shear  strength  properties  attributed  to  the  four  layers,  and 
utilized  in  calculations,  account  for  the  stair-step  deceleration  his- 
tory calculated  as  the  projectile  penetrates  layers  3 and  4.  Measured 
deceleration  levels  in  layers  1 and  2,  as  indicated  by  the  data  trace, 

the  smoothed  (by  a = f ^ a(t)dt / j'  dt)  data  curve,  and  the  effective  con- 

0 0 

stant  deceleration  over  each  designated  layer  are  greater  than  calculated. 
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Comparison  of  Predicted  and  Measured  Projectile  Deceleration  versus  Depth. 
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Comparing  calculated  and  measured  average  deceleration  levels 
(Table  V-l)  for  each  of  the  first  3 layers  it  is  seen  that  predicted 
values  are  23 i 31  and  13$  lower,  respectively  than  observed.  For  layer 
4 uhe  difference  is  considerably  greater  and  the  reasons  for  this  are 
discussed  further  below. 


5.  * 


Since  the  average  decelerations,  predicted  by  the  no  friction  cal- 
culation, are  less  than  those  observed  one  may  speculate  whether  or  not 
''  the  differences  are  a measure  of  the  friction  forces  operative.  For 

layers  1,  2 and  3)  it  is  seen  (Table  V-l)  that  the  difference  is  greatest 
for  ••■he  layer  2 of  highest  shear  strength  and  least  for  the  layer  3 of 
■ lowest  shear  strength.  It  is  apparent  that  the  differences  between 

j.  experimental  and  predicted  average  decelerations  are  in  monotonic  cor- 

* respondence  to  shear  strengths  of  the  materials  at  the  normal  stress 

1 

j experienced  on  the  projectile  nose  (Figures  A-98,  A-ll6  and  A-134) . 

I 

i This  suggests  that  a form  of  the  friction  law  involving  the  second 

i 

1 invariant  of  the  deviator  stresses,  such  as  t»  a t (t  * frictional 

l 1 1 

1 

3 stress  along  soil  projectile  interface,  t - octahedral  shear  stress), 

< is  not  inconsistent  with  data.  Similarly,  the  data  of  Table  V-l  sug- 

i gests  that  a friction  law  related  to  projectile  velocity  is  probably  not 

.ji  correct. 

/i 

j;.  It  is  of  interest  to  explore  the  possible  sources  of  discrepancy 

p 

? between  calculations  and  experiment  to  determine  whether  or  not  more 

h 

ji  definitive  information  regarding  friction  forces  may  be  obtained 
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Table  V-l 


Difference  Between  Predicted  and  Measured  Average  Projectile 
Decelerations  in  the  Four  Layers  Assumed  for  the 
Watching  Hill  Medium 


Average  Deceleration  (g's) 
Calculated*  Observed 


Difference 


50 

65 

23 

55 

80 

31 

35 

4o 

13 

< 20 

> 60 

> 67 

Layer 

1 

2 

3 

4 

* No  Friction 

in  order  to  refine  or  perfect  the  prediction  technique.  Sources  of 
discrepancy  include  errors  in  the  deceleration  measurement,  errors  in 
laboratory  measurements  of  material  properties,  approximations  or  assum- 
ptions made  in  calculation  of  the  penetration  problem  and  possible  in- 
accuracies associated  with  the  computational  method.  Approximations  and 
assumptions  made  in  calculation  of  the  penetration  problem  include  neg- 
lect of  variations  in  lithostatic  pressure  with  depth,  neglect  of  friction 
forces,  layering  in  target  and  use  of  the  soil  cap  material  response  model 
to  represent  the  target  medium  Inaccuracies  in  the  computational  method 
involve  size  of  mesh  and  the  closeness  with  which  the  cap  model  is  made 
to  represent  laboratory  data  on  material  properties  of  the  medium. 

While  the  uncertainty  associated  with  each  of  these  sources  of 
error  is  not  known,  it  is  expected  that  taken  in  combination  ail  the 
sources  of  error  from  both  experiment  and  calculation  are  great  enough 
to  account  for  the  13  to  31 % discrepancies  noted  between  predictions  and 
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observations.  As  a result,  no  definitive  conclusions  can  be  drawn  re- 
garding the  absolute  magnitude  of  friction  forces  contributing  to  de- 
celeration of  the  projectile  in  layers  1,  2 and  3*  However,  it  appears 
that  the  friction  forces  are  most  likely  less  than  50%  of  the  total 

axial  force  acting  on  the  projectile.  California  Research  and  Technology 
20 

predictions,  obtained  with  friction  on  the  projectile  specified  by 

T 


= 0.6  \ jl  , yielded  deceleration  amplitudes  of  110  g and  125  g for 
x max 

layers  1 and  2 respectively.  Here,  yjjl  is  the  function  describing 
the  failure  envelope,  evaluated  at  the  normal  stress  acting  on  the  pro- 
jectile. The  contribution  from  friction  to  the  total  axial  retarding 
force  was  about  55%  and  56%,  respectively.  This  amount  of  friction  is 
clearly  too  great-,,  producing  decelerations  larger  than  observed.  At 
most,  the  amount  of  friction  for  layers  1,  2 and  3 is  that  suggested 
by  differences  indicated  in  Table  V-l. 

The  difference  between  measured  and  calculated  decelerations  in 
layer  4 material  is  likely  greater  than  can  be  accounted  for  by  uncer- 
tainties in  the  experimental,  data  on  deceleration  and  in  the  computation- 
al method.  The  large  discrepancy  is  felt  to  be  a consequence  of  material 
property  data  and  the  resultant  material  response  model  which  is  not 

representative  of  the  actual  material  of  layer  4.  Credence  to  this 

5 

supposition  is  given  by  the  full  experimental  deceleration  record 
shown  in  Figure  V-3«  At  a depth  of  7.2  m,  where  layer  4 is  assumed  to 
begin,  projectile  deceleration  begins  to  increase  and  from  7.2  m to  10 
m the  projectile  experiences  an  increase  in  retarding  force  which  implies 
a layer  of  greater  shear  strength  than  that  of  layer  3*  Below  a depth 
of  10  m the  projectile  enters  another  layer  of  lesser  strength  than  the 
last  (but  still  greater  than  layer  3 strength)  and  maintains  a nearly 
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constant  deceleration  of  slightly  more  than  60  g's.  These  observations 
indicate  that  the  target  medium  would  have  been  modeled  better  with  five 
layers  rather  than  four  and  also  that  the  material  properties  used  for  layer 
four  are  not  representative  of  any  of  the  layers  assumed  for  the  penetrated 
area  of  the  Watching  Hill  medium.  In  retrospect,  it  is  apparent  that  a 
prior  choice  of  four  layers  for  use  in  calculations  was  not  the  best. 
Considerably  improved  results  for  predicted  deceleration  history  would 
be  obtained  in  a recalculation  of  the  problem  by  using  five  layers  rather 
than  four,  by  taking  layer  four  to  be  much  stiff er  in  shear  than  that 
used  in  the  first  calculation,  and  also  by  employing  a small  amount  of 
friction. 

Further  comparisons  of  calculated  projectile  motion  with  information 
derived  from  deceleration  data  are  shown  in  Figures  V-4  through  V-8. 

Of  particular  interest  among  these  is  Figure  V-7  for  deceleration  versus 
velocity.  Within  a given  layer  the  calculated  curve  is  independent  of 
projectile  velocity  indicating  a constant  retarding  force  within  each 
uniform  layer  of  material.  This  result  for  the  Watching  Hill  medium  seems 
to  confirm  the  force  law  relationship  proposed  by  Robins  and  Euler\  as 
opposed  to  those,  such  as  Resal's,^  which  require  an  explicit  dependence 
on  projectile  velocity. 

Target  Medi un  Respons e_ 

To  measure  medium  response  during  projectile  penetration,  six  lithium 
niobate  stress  gauges  were  located^  at  positions  as  indicated  in  Figure 
V-9a.  Figure  V-9b  illustrates  the  orientation  of  four  penetration  experi- 
ments about  the  stress  transducer  locations  and  Table  V-3  lists  particu- 

5 

lars  for  those  four  penetrator  experiments. 
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Figure  V-5 

Comparison  of  Predicted  and  Measured  Projectile  Velocity  versus  Depth. 
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Figure  V-8 

Projectile  Kinetic  Energy  versus  Time,  Measured  and  Predicted. 
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Figure  V-9a 

Vertical  Location  of  Earth  Stress  Transducers  and  Orientation  of 
Stress  Gauge  Holes  Relative  to  Penetrator  Impact  Positions. 
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Table  V-3 

Characteristics  of  Projectiles  Penetrating 
Medium  in  Neighborhood  of  Stress  Gauges 


Test  No. 
1 
2 

3 

4 


Projectile 

Mass  (Kg)  Diameter  (m) 


127 

181 

12? 

127 


0.152 

0.165 

0.152 

0.152 


Velocity  (m/s) 
174 
125 
122 
156 


tv 


No  detectable  signals  were  sensed  by  any  stress  transducers  which 
were  more  than  1.2  m horizontal  distance  from  the  projectile  path. 
Excellent  stress  histories  were  measured  for  penetrator  test  number  4. 
Also,  the  gauge  in  hole  B measured  a signal  from  penetrator  test  number 
2.  Data  were  obtained  from  penetrator  test  number  3 also,  but  were  un- 
fortunately lost  during  playback  of  recording  tapes. 

Measured  stress  histories  are  compared  with  those  calculated  in 
Figures  V-10  through  V-15,  where  it  is  seen  that  amplitudes,  pulse 
shapes  and  arrival  times  are  all  in  generally  good  agreement  with  predic- 
tions. Zero  time  in  Figures  V-10  through  V-15  is  when  the  projectile 
is  0.3  m above  target  surface."^  Peak  values  of  stresses  are  summarized 
in  Table  V-4  where  uncertainties  in  measured  peak  stresses  associated 
with  base  line  shift  are  also  tabulated.  Error  in  peak  stresses  is  11% 
or  less  for  stress  transducers  closest  to  the  impact  axis  (0.45  m) . At 
greater  distances  from  the  impact  axis,  error  due  to  base  line  shift 
increases  appreciably  as  the  magnitude  of  the  signal  goes  down.  Errors 
from  other  possible  sources  are  difficult  to  estimate.  The  extent  to 
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which  earlier  penetration  experiments  disturb  the  medium  and  possibly 
influence  the  accuracy  of  stress  history  measurements  on  subsequent 
experiments  cannot  be  ascertained  with  available  data.  The  disturbance 

I 

produced  by  penetrator  test  number  3 (Figure  V-9b)  in  terms  of  calculated 
peak  displacement  at  the  gauge  positions  amounts  to  about  10  mm.  Further 

» 

experiments  are  necessary  to  evaluate  whether  or  not  this  magnitude  dis- 
turbance affects  the  measurements  made  on  penetrator  test  number  4.  It 
is  clear  that  penetrator  tests  numbers  1 and  2 produced  relatively  mucn 
smaller  disturbances  at  the  gauges. 

Shown  in  Figure  V-l6  are  peak  pressure  and  radial  stress  data  from 
gauges  in  hole  A plotted  as  a function  of  gauge  depth.  Also  shown  are 
the  calculated  curves  for  peak  pressure  and  radial  stress.  Except  for 
> the  datum  from  gauge  A-8,  all  the  data  fall  within  a band  bounded  by  the 

calculated  curves.  Also,  all  data  are  closer  to  the  predicted  mean 
stress  curve  than  to  the  radial  stress  curve.  The  datum  at  3*76  m from 
the  radial  stress  gauge  (A-5)  appears  to  be  essentially  identical  to  the 
pressure  measurement  from  gauge  A-7  suggesting  that  both  gauges  are 
measuring  the  same  quantity.  Lithium  niobate  transducers  for  measure- 
ment of  total  stress  components  are  still  under  development.  More  work 
is  necessary  to  establish  their  capability  in  this  mode.  Mean  stress 

q 

gauges  of  lithium  niobate,  however,  have  received  much  more  attention 
2] 

and  use  and  should  be  considered  the  more  reliable  transducer  in  examin- 
ation of  data  from  the  penetration  experiments , 

‘ Attenuation  of  peak  pressure  and  peak  radial  stress  with  radial 

position  from  impact  axis  is  illustrated  in  B'igure  V-17-  Predicted 
. curves  shown  are  for  stress  attenuation  at  a depth  of  3*6  m and  may  be 

compared  with  data  at  that  depth  (circle  and  triangles).  It  is  seen 
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Figure  V-16 

Predicted  Peak  Radial  and  Peak  Mean  Stress  versus  Depth  at 
Constant  Radial  Position  R = 0.45  Compared  with  Measurements. 
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that  the  3*6  m data  compare  reasonably  well  with  the  predicted  mean  stress 

attenuation  curve  and  that  peak  stress  amplitude  falls  off  rapidly  (faster 

-2 

than  r ) with  distance  away  from  the  projectile  path  axis. 

In  general,  it  is  felt  that  the  soil  stress  measurement  technique 
has  been  demonstrated  and  that  the  information  obtained  provides  con- 
fidence regarding  validity  of  the  material  property  data  and  material 
response  models  used  in  calculations  for  layers  1,  2 and  3«  The  satisfac- 
tory agreement  between  measurement  and  calculation  of  st-ess  pulse  shape 
and  amplitude  also  verifies  methods  employed  in  the  computational  tech- 
nique. Without  the  earth  stress  data  much  less  credibility  might  have 
been  attached  to  results  of  computations.  Penetration  tests  (numbers  2 
and  4) , which  provided  the  earth  stress  history  data  were  conducted  with 
projectiles  of  12?  kg  mass  and  diameter  0.152  m rather  than  the  l8l  kg, 
O.165  m diameter  projectiles  assumed  in  calculations.  While  another 
calculation  would  be  required  for  the  slightly  smaller  projectile  in 
order  to  definitely  establish  differences,  it  is  felt  that  stress  history 
results  for  the  smaller  projectile  would  be  similar  to  those  of  the 
larger  projectile  within  the  uncertainties  that  exist. 
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SECTION  VI 


SUMMARY  AND  CONCLUSIONS 

All  phenomena  connected  with  producing  significant  loads  on  the 
projectile  occur’  in  a relatively  small  region  of  the  target  medium  about 
the  projectile  nose.  Shown  in  Figure  VI-1  are  three-dimensional  plots 
of  radial  and  vertical  soil  velocity  as  a function  of  depth  in  target 
and  distance  from  projectile  path  axis  at  a time  when  the  projectile 
nose  is  at  4.4  m depth  (see  also  Figures  A-122  and  A-123) . It  is  seen 
that  the  region  of  greatest  disturbance  in  the  soil  is  small  and  confined 
to  the  immediate  vicinity  of  the  nose.  Figure  VI-2a  also  illustrates 
how  rapidly  the  disturbance  in  soil  vanishes  with  distance  from  the 
projectile.  With  the  projectile  nose  tip  at  4.4  m,  at  a depth  of  4.2  m 
(0.2  in  behind  projectile  nose)  soil  pressure  is  about  3 MPa  (30  bars, 

450  psi)  at  a distance  of  slightly  less  than  one  projectile  radius  from 
the  projectile  axis.  At  a distance  2.5  projectile  diameters  from  the 
axis,  the  soil  pressure  is  down  an  order  of  magnitude  to  0.3  MPa  and  at 
5 projectile  diameters  another  order  pf  magnitude  reduction  is  found. 
Figure  VI-2b  indicates  in  yet  another  way  the  highly  isolated  nature  of 
the  projectile  generated  disturbance  in  the  soil.  For  the  zones  of  the 
computational  mesh  originally  on  the  symmetry  axis  and  later  adjacent 
to  the  penetrator,  mean  stress  is  plotted  as  a function  of  vertical 
position  relative  to  nose  tip  location.  The  pressure  profile  ahead  of 
and  along  the  projectile  surface  is  shown  in  Figure  VI-2b  for  the  penetra- 
tor in  layer  1 (nose  tip  at  2 m depth)  and  in  layer  2 (4.4  m depth). 

For  materials  similar  to  those  of  the  Watching  Hill  target  medium,  it  is 
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seen  that  essentially  al],  of  the  disturbance  produced  in  the  soil  occurs 
within  a cylindrical  volume  with  radius  equal  to  3 projectile  diameters. 
For  purposes  of  calculating  loads  on  penetrators , radial  extent  of  the 
target  need  not  be  more  than  about  5 penetrator  diameters  in  materials 
like  those  of  this  calculation. 

Figure  VI-3  shows  the  path  a soil  particle  follows  in  the  stress 
plane,  t vs.  p.  The  particle  of  Figure  VI-3  is  identified  by  its 
initial  coordinates  Z = 3-6  m,  R = 0.15  m and  so  is  located  in  layer  2 
material  at  a distance  less  than  one  projectile  diameter  from  the  impact 
axis.  During  loading  the  particle  follows  the  path  shown.  On  unloading 
the  particle  finds  itself  very  close  to  the  failure  envelope  and  follows 
it  to  zero  pressure.  A close  examination  of  the  calculation  reveals  that 
the  stress  state  for  the  particle  should  actually  be  at  the  intersection 
of  the  cap  and  the  failure  envelope:  the  fact  that  it  does  not  quite 

coincide  with  this  point  is  a consequence  of  a slight  "over-shrinking" 
of  the  cap  caused  by  the  finite  difference  approximation  used  for  the 
constitutive  .equations.  Figures  VI-2b  and  VI-3  taken  together  indicate 
that  accurate  modeling  of  the  failure  envelope  and  unloading  behavior 
for  the  material  are  important  parts  of  the  calculation.  Similarly, 
Figure  VI-4  illustrates  the  path  of  the  same  soil  particle  in  the  strain 
invariant  plane,  C vs.  T],  where  I g is  the  second  moment  of  the  strain 
deviators  and  T|  is  volumetric  strain.  After  unloading,  the  particle 
finds  itself  with  a large  deviatoric  strain  and  at  a density  less  than 
its  initial  density.  These  paths  experienced  by  soil  particles  as  a 
result  of  target  penetration  by  projectiles  are  considerably  different 
from  paths  attained  in  conventional  uniaxial  strain  and  triaxial  stre  j 
experiments . Strain  energy  density  of  a particle  within  the  region  of 
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soil  highly  disturbed  by  the  projectile  is  of  the  order  of  0.1  to  0.5  MJ/ 

for  the  Watching  Hill  medium.  Consequences  of  the  dilatancy  inherent 

in  the  soil  cap  model,  as  indicated  in  Figure  VI-4,  appear  to  be  manifested 

in  two  ways.  These  are  lack  of  a definitive  soil  separation  point  on  the 

projectile  and  the  size  of  the  projectile-produced  borehole.  Results  of 

the  rigid  body  calculation  indicate  that  the  diameter  of  the  borehole  is 

never  greater  than  the  projectile  diameter  and  may  in  fact  be  less. 

Accurate  determination  of  borehole  diameter  behind  the  projectile  is  not 

possible  owing  to  the  continual  rezoning  to  a coarser  mesh  in  this  region. 

It  is  known,  however,  that  soil  rebound  occurs  and  boreholes  have  been 
22 

observed  whose  final  diameters  are  less  than  the  projectile  diameter. 
Unfortunately,  no  attempt  was  made  to  measure  borehole  sizes  for  the 
penetrator  experiments  at  Watching  Hill.  The  rigid  body  calculations 
also  suggest  that  the  Watching  Hill  soil  does  not  separate  from  the 
projectile  but  appears  to  remain  attached  throughout  the  entire  length 
of  the  projectile.  This  occurs  even  though  the  normal  stress  on  the 
projectile  is  significantly  larger  than  zero  only  on  the  nose.  (What 
appears  to  be  a soil  separation  position  in  Figures  A-97>  A-lOo,  A-115, 
A-124  and  A-133  is  interpreted,  rather,  to  be  a perturbation  in  the  cal- 
culation resulting  from  the  transition  from  a region  of  fine  to  coarser 

zoning) . That  soil  separation  or  detachment  from  the  projectile  may  not 

22 

occur  has  been  observed  experimentally.  Those  experimental  results, 
however,  cannot  distinguish  between  soil  detachment  on  the  projectile 
nose  with  later  reattachment  and  complete  lack  of  soil  separation. 

To  investigate  this  point  further,  additional  calculations  should 
be  performed  which  include  friction  and  a response  model  such  as  the 
cap  model  which  demonstrates  dilatancy.  As  appears  to  be  the  case  here 
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in  calculations  with  no  friction,  dilatancy  in  the  response  model  pre- 
vents or  retards  soil  detachment  from  the  projectile  surface.  By  the 
inclusion  of  friction  in  such  calculations  even  a small  frictional  stress, 
acting  over  the  large  area  of  the  projectile  afterbody,  could  contribute 
a large  axial  retarding  force  comparable  to  or  possibly  greater  than 
that  contributed  by  the  nose  area  of  the  projectile.  Results  from  such 
a calculation,  when  compared  with  projectile  deceleration  data,  would 
imply  a smaller  friction  coefficient  than  one  inferred  from  calculations 
with  response  models  having  no  dilatancy  in  which  soil  separation  might 
occur,  for  example,  at  the  position  on  the  projectile  where  the  nose 
meets  the  afterbody. 

Figure  VI-5  illustrates  the  partitioning  of  energy  to  the  penetrated 
medium  at  the  expense  of  projectile  kinetic  energy.  From  the  zero  fric- 
tion calculation,  it  is  apparent  that  projectile  kinetic  energy  losses 
manifest  themselves  primarily  in  the  form  of  plastic  work  done  in  defor- 
mation of  the  penetrated  medium.  Energies  associated  with  motion  of  the 
soil  and  with  recoverable  elastic  work  are  very  small,  each  being  on  the 
order  of  'l/lO  or  less  of  the  plastic  work  (internal  energy  = plastic  work  + 
elastic  work) . 

The  calculated  deceleration  history  indicates  that  sudden  changes 
in  the  problem  occur  only  when  the  nose  region  of  the  projectile  is  entering 
a new  layer  of  material.  Also,  Figure  VI-5  shows  that  the  dominant  mode 
of  the  rate  of  energy  transfer  from  projectile  to  soil  is  nearly  constant 
through  a layer  in  which  there  is  little  velocity  change.  Moving  pictures 
made  of  numerous  plots  of  normal  stress  distribution  and  three-dimensional 
plots  of  pressure  vs.  axial  and  radial  position,  for  example,  confirm 
that  the  process  being  calculated  is  primarily  quasi-static.  Thus,  the 
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method,  of  "skipping"  through  layers,  and  calculating  only  transition 

20 

periods,  as  employed  by  California  'Research  and  Technology,  appears 
to  be  a useful  means  of  performing  deep  penetration  calculations  with 
large-scale  wave  propagation  codes  with  a reasonable  degree  of  economy 
although  it  may  be  necessary  to  perform  multiple  calculations  in  layers 
that  produce  large  velocity  decreases.  It  also  appears,  however,  that 
detailed  analysis  of  the  results  of  this  and  similar  calculations  should 
lead  to  the  development  of  much  less  sophisticated  methods  which  would 
also  yield  necessary  information  for  depth  prediction  and  projectile  design. 

During  the  deformable  projectile  calculation,  the  maximum  radial 
displacement  that  was  observed  for  the  projectile  surface  was  30 
Although  this  displacement  was  too  small  to  alter  the  effective  shape 
of  the  vehicle,  fluctuations  in  displacement  were  sufficient  to  produce 
large  charges  in  normal  stress  at  the  vehicle-soil  interface,  as  shown 
in  Figures  B-3?  through  B-39*  Thus,  the  dominant  effect  of  the  rigid- 
body  approximation  was  to  eliminate  transient  variations  in  the  normal 
stress. 

The  agreement  in  predicted  soil  histories  between  the  two  calcula- 
tions demonstrates  that,  for  the  particular  material  model  used,  these 
transients  are  damped  out  less  than  20  mm  from  the  vehicle  and  do  not 
affect  the  predicted  wave  motion  in  the  soil.  However,  the  deformable 
body  calculation  predicts  a deceleration  about  forty  percent  larger  than 
that  calculated  for  a rigid  body.  This  comparison  demonstrates  that 
the  deceleration  is  determined  by  the  soil  response  in  regions  much 
closer  to  the  vehicle  than  the  locations  at  which  stress  measurements  were 
made.  The  motion  is  similar  to  aerodynamic  behavior,  with  the  drag  de- 
termined by  the  response  in  a narrow  region,  or  boundary  layer. 


75 


I 


sihlil: iiAff&lAS-yf  A3 


ifll  livS  L-  ■ 


/’ 

/ 

/ 

/. 

The  differences  between  the  results  of  the  two  calculations  are 
related  to  their  description  of  this  boundary  layer.  The  inclusion  of  the 
45°  nose  tip  results  in  a stress  increase  of  less  than  5 MPa,  effective  over 

i -4  2 

a cross-sectional  area  of  about  4 x 10  m , accounting  for  a deceleration 
difference  of  about  1 g.  Of  greater  significance  is  the  procedure  for 
obtaining  the  stress  at  the  interface.  In  the  calculations,  stress  his- 
tories are  recorded  only  at  zone  centers  and,  therefore,  not  at  the 
interface.  For  the  rigid  body  calculation,  the  stress  in  the  nearest 
soil  zone  was  used  to  determine  the  decelerating  force,  as  discussed 
in  Section  III.  In  the  deformable  body  calculation,  no  interface  stress 
or  force  was  explicitly  obtained;  the  sliding  interface  procedure  results 
in  an  effective  normal  stress  that  is  between  the  stress  in  the  soil  zone 
and  the  (usually  higher)  stress  in  the  nearest  vehicle  zone. 

While  many  questions  remain  to  be  investigated,  this  first  attempt 
to  predict  projectile  motion  and  target  material  response  has  been 
highly  successful  and  most  encouraging.  Given  data  describing  properties 
of  each  layer  in  the  target  medium,  realistic  soil  response  models  were 
constructed  and  utilized  in  the  finite  difference  code  computations. 

Resulting  predictions  of  projectile  deceleration  and  of  stress  history 
at  several  points  in  the  soil  generally  agreed  well  with  measurements, 
lending  credence  to  the  method  of  analysis.  Consequently  greater  con- 
fidence may  now  be  placed  in  prediction  of  nonmeasurable  quantities, 
such  as  normal  stress  distribution  on  projectile  surface,  which  are 
necessary  for  design  of  earth  penetrating  projectiles. 
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List  of  Illustrations  Appendix  A 


Figure  A-l  Projectile  Deceleration  vs  Time. 

Figure  A -2  Projectile  Velocity  vs  Time. 

Figure  A -3  Projectile  Depth  vs  Time. 

Figure  A-k  Projectile  Deceleration  vs  Depth. 

Figure  A-5  Projectile  Velocity  vs  Depth. 

Figure  A-6  Projectile  Kinetic  Energy/Radian  vs  Time, 

Figure  A -7  Soil  Kinetic  Energy/Radian,  Due  to  Radial  Velocity,  vs  Time. 

Figure  A-8  Soil  Kinetic  Energy/ Radian,  Due  to  Axial  Velocity,  vs  Time. 

Figure  A -9  Total  Soil  Kinetic  Energy/Radian  vs  Time. 

Figure  A -10  Legend  for  Point  History  Plots  from  Rigid  Projectile  Calculation. 

Figure  A -11  Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  3> 

Z = 3.6  m,  E = 0.15  m, 

Figure  A -12  Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  3> 

Z = 3.6  m,  R = 0.15  m. 

Figure  A -13  Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  3> 

Z = 3.6  m,  R = 0.15  m. 

Figure  A-lk  Axial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  3> 

Z = 3.6  m,  R = 0.15  m. 

Figure  A -15  Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  3,  Z = 3.6  m, 
R = 0.15  m. 

Figure  A-l6  Deviator  Stress  , J2,  History  (Pa)  for  Particle  at  Point  3>  Z = 
3.6  m,  R = 0.15  m (Octahedral  Shear  Stress  = \/2/3(j2) ) . 

Figure  A -17  Radial  Stress  (Pa)  History  for  Particle  at  Point  3>  Z = 3.6  m, 

R = 0.15  m. 

Figure  A-18  Axial  Stress  (Pa)  History  for  Particle  at  Point  3>  Z - 3.6  m, 

R “ 0.15  m. 
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Figure  A-ig 


Figure  A -20 


Figure  A -21 


Figure  A -22 


Figure  A -23 


Figure  A -24 


Figure  A-25 


Figure  A-26 


Figure  A -27 


Figure  A -28 


Figure  A -29 


Figure  A -30 


Figure  A -31 


Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  4, 

Z = 0.4  in,  R = 0.15  tn. 

Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  4, 

Z = 0.4  m,  R = 0.15  m. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  4, 

Z = 0.4  m,  R = 0.15  m* 

Axial  Displacement  (m)  versus  Time  (p)  for  Particle  at  Point  4, 

Z = 0.4  m,  R = 0.15  m- 

Pressure  (pa)  versus  Time  (s)  for  Particle  at  Point  4,  Z = 0.4  rn, 
R - 0.15  m. 

Deviator  Stress,  J2,  History  (Pa)  for  Particle  at  Point  4, 

Z = 0.4  m,  R - 0.15  m (Octahedral  Shear  Stress  - V^/3(J2) ) . 
Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  5> 

Z = 0.1  m,  R = 0.15  ni. 

Axial  Velocity  (m/s)  versus  Time  (s)  for  particle  at  Point  5> 

Z = 0.1  m,  R = 0.15  m. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  5> 

Z = 0.1  m,  R = 0.15  m. 

Axial  Displacement  (in)  versus  Time  (s)  for  Particle  at  Point  5> 

Z = 0.1  m,  R = 0.15  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  5>  Z,  = 0.1  m, 
R = 0.15  m. 

Deviator  Stress,  J2,  History  (Pa)  for  Particle  at  Point  5> 

Z = 0.1  m,  R = 0,15  m (Octahedral  Shear  Stress  =V2/3(J2)V. 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  7> 

Z = 5.4  m,  R = 0.45  m. 
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Figure  A-32 
Figure  A -33 
Figure  A-34 
Figure  A-35 
Figure  A -36 
Figure  A -37 
Figure  A -38 
Figure  A-39 
Figure  A -1+0 
Figure  A-4l 
Figure  A -42 
Figure  A-43 

Figure  A -44 


Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  7, 

Z = 5.4  m,  R - 0.45  m. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  7, 

Z = 5.4  m,  R = 0,45  m. 

Axial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  7 > 

Z = 5.4  m,  R - 0.45  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  7>  Z,  = 5.4  m, 
R = 0.4-5  m. 

Deviator  Stress,  J2,  History  (Pa)  for  Particle  at  Point  7, 

Z = 5.4  m,  R = 0.45  m (Octahedral  Shear  Stress  - \/2/3(J2) ) . 
Radial  Stress  (Pa)  History  for  Particle  at  Point  7>  Z = 5*4  m, 

R = 0.4?  m. 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  8, 

Z = 3.6  m,  R = 0.45  m. 

Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  8, 

Z = 3.8  m,  R = 0.45  m. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  Q, 

Z = 3.6  m,  R = 0.45  m. 

Axial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  8, 

Z = 3.8  m,  R = 0.45  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  8,  Z = 3.6  m, 
R = 0.45  m. 

Deviator  Stress,  J2,  History  (Pa.)  for  Particle  at  Point  8, 

Z = 3*6  m,  R = 0.45  m (Octahedral  Shear  Stress  = v£73(j2)). 
Radial  Stress  (Pa)  History  for  Particle  at  Point  8,  Z = 3.6  in, 

R = 0.45  m. 
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Figure  A -45 


Figure  A -46 


Figure  A -47 


Figure  A -48 


Figure  A -49 


Figure  A-50 


Figure  A-51 


Figure  A -52 


Figure  A -53 


Figure  A -54 


Figure  A-55 


Figure  A-56 


Figure  A -57 


Axial  Stress  (Pa)  History  for  Particle  at  Point  8,  Z = 3.6  in, 

R = 0.4-5  m. 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  9, 

Z = 1.8  m,  R - 0.45  m. 

Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  9, 

Z 1.8  m,  R = 0.45  m. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  9, 
Z = 1.8  m,  R = 0.4-5  m. 

Axial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  9, 

Z = 1.8  in,  R = 0.45  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  9,  Z = 1.8  m, 
R = 0.4-5  m. 

Deviator  Stress,  J2,  History  (Pa.)  for  Particle  at  Point  9, 

Z = 1.8  m,  R = 0.45  m (Octahedral  Shear  Stress  =V2/3(J2) ) . 

Radial  Stress  (Pa)  History  for  Particle  at  Point  9,  Z = 1.8  m, 

R » 0.45  m. 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  10, 

Z = 0.4  m,  R a 0.45  m. 

Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  10, 

Z - 0.4  m,  R = 0.45  m. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  10, 
Z = 0.4  m,  R = 0.45  m. 

Axial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  10, 

Z = 0.4  m,  R = 0.45  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  10, 

Z = 0.4  m,  R = 0.45  m. 
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Figure  A -58 
Figure  A-59 
Figure  A-60 
Figure  A-6l 
Figure  A -62 
Figure  A -63 
Figure  A -64- 
Figure  A-65 
Figure  A -66 
Figure  A-67 
Figure  A -68 
Figure  A -69 

Figure  A-70 


Deviatur  Stress,  J2,  History  (Pa)  for  Particle  at  Poirr  " . Z - 
0.4  m,  R = 0.45  m (Octahedral  Shear  Stress  = >/2/3(J2)). 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  11, 

Z = 0.1  m,  R = 0.45  m. 

Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  11, 

Z = 0.1  m,  R = 0.45  m. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  11, 
Z = 0.1  m,  R - 0.4-5  m. 

Axial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  11, 

Z = 0,1  m,  R - 0.45  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  11,  Z = 

0.1  m,  R = 0.45  m. 

Deviator  Stress,  J2,  History  (Pa)  for  Particle  at  Point  11,  Z = 
0.1  m,  R = 0.45  m (Octahedral  Shear  Stress  = V2/3(j2) ) . 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  12, 

Z = 3.6  m,  R = 0.9  m. 

Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  12, 

Z = 3«6m,  E = 0,9  11. 

Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  12, 
Z = 3.6  m,  R = 0.9  m. 

Axial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  12, 

Z = 3.6  m,  R = 0.9  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  12,  Z = 

3.6  ra,  R = 0.9  m. 

Deviator  Stress,  J2,  History  (Pa)  for  Particle  at  Feint  12,  Z = 
3.6  m,  R = 0.9  m (Octahedral  Shear  Stress  = \f2j^{S2) ) . 
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Figure  A-71 

Figure  A -72 

Figure  A -73 

Figure  A-74 

Figure  A -75 

Figure  A-76 

Figure  A -77 

Figure  A -78 

Figure  A-79 

Figure  A-80 

Figure  A-8l 

Figure  A-82 
Figure  A-83 
Figure  A-84 


Radial  Stress  (Pa)  History  for  Particle  at  Point  12,  Z = 3.6  m, 

R = 0.9  m. 

Axial  Stress  (Pa)  History  r>r  Particle  at  Point  12,  Z = 3.6  m, 

R = 0.9  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  13, 

Z = 5.4  m,  R = 1.5  m. 

Radial  Stress  (pa)  History  for  Particle  at  Point  13,  Z = 5.4  m, 

R = 1.5  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  14,  Z = 3.6  m, 
R = 1.5  m. 

Radial  Stress  (Pa)  History  for  Particle  at  Point  14,  Z = 3.6  m, 

R = 1.5  m. 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  15,  Z = 1.8  m, 
R = 1.5  m. 

Radial  Stress  (Pa)  History  for  Particle  at  Point  15,  Z = 1.8  m, 

R = 1.5  m. 

Lagrangian  coordinates  when  projectile  has  penetrated  to  depth 
of  0.1  m. 

Normal  stress  (Pa)  along  surface  of  projectile  when  projectile 
is  at  0.1  m depth. 

Tangential  stress  (Pa)  along  surface  of  projectile  when  projectile 
is  at  0.1  m depth. 

Contours  of  Axial  Stress  in  Target.  Projectile  at  0.1  in  Depth. 
Contours  of  Radial  Stress  in  Target.  Projectile  at  0.1  m Depth. 
Contours  of  Hoop  Stress  in  Target.  Piujectile  at  0.1  m Depth. 
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Figure  A-85  Contours  of  Radial-Vertical  Shear  Stress  in  Target.  Projectile 
at  0.1  m Depth 

Figure  A-86  Contours  of  Vertical  Velocity  in  Target.  Projectile  at  0.1  m Depth. 

Figure  A-87  Contours  of  Radial  Velocity  in  Target.  Projectile  at  0.1  m Depth. 

Figure  A-83  Lagrangian  Coordinates  When  Projectile  has  Penetrated  to  Depth 
of  0.4  m. 

Figure  A-89  Normal  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  0.4  m Depth. 

Figure  A-90  Tangential  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  0.4  m Depth. 

Figure  A-91  Contours  of  Axial  Stress  in  Target.  Projectile  at  0.4  m Depth. 

Figure  A-92  Contours  of  Radial  Stress  in  Target.  Projectile  at  0.4  m Depth. 

Figure  A-93  Contour’s  of  Hoop  Stesss  in  Target.  Projectile  at  0.4  m Depth. 

Figure  A-94  Contours  of  Radial-Vertical  Shear  Stress  in  Target.  Projectile 
at  0.4  m Depth. 

Figure  A-95  Contours  of  Vertical  Velocity  in  Target.  Projectile  at  0.4  m Depth. 

Figure  A-96  Contours  of  Radial  Velocity  in  Target.  Projectile  at  0.4  m Depth. 

Figure  A-97  Lagrangian  Coordinates  When  Projectile  has  Penetrated  to  Depth  of 
2.0  m.  Radial  Dimensions  Doubled  for  Clarity. 

Figure  A-98  Normal  Stress  (Pa)  along  Surface  of  Projectile  When  Projectile 
is  at  2.0  ra  Depth. 

Figure  A-99  Tangential  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  2.0  m Depth. 

Figure  A-100  Contours  of  Axial  Stress  in  Target.  Projectile  at  2.0  m Depth. 

Figure  A-101  Contours  of  Radial  Stress  in  Target.  Projectile  at  2.0  m Depth. 

Figure  A-102  Contours  of  Hoop  Stress  in  Target.  Projectile  at  2.0  m Depth. 
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Figure  A- 103  Contours  of  Radial  Vertical  Shear  Stress  in  Target.  Projectile 
at  2.0  m Depth, 

Figure  A-104  Contours  of  Vertical  Velocity  in  Target.  Projectile  at  2.0  ra  Depth. 

Figure  A-105  Contours  of  Radial  Velocity  in  Target.  Projectile  at  2.0  m Depth. 

Figure  A-lCo  Lagrangian  Coordinates  When  Projectile  has  Penetrated  to  Depth 
of  2.8  in.  Radial  Dimensions  Doubled  for  Clarity. 

Figure  A -107  Normal  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  2.8  m Depth. 

Figure  A-108  Tangential  Stress  (Pa)  Along  Surface  of  Projectile  when  Projectile 


Figure  A-I09 
Figure  A-110 
Figure  A-lll 


Figure  A -112 


is  t 2.8  m Depth. 

Jontours  of  Axial  Stress  in  Target.  Projectile  at  2.8  m Depth. 
Contours  of  Radial  Stress  in  Target.  Projectile  at  2.8  m Depth, 
Contours  of  Hoop  Stress  in  Target.  Projectile  at  2.8  m Depth. 
Contours  of  Radial  Vertical  Shear  Stress  in  Target.  Projectile 
at  2.8  m Depth. 


Figure  A-113  Contours  of  Vertical  Velocity  in  Target.  Projectile  at  2.8  m Depth. 
Figure  A-114  Contours  of  Radial  Velocity  in  Target.  Projectile  at  2.8  m Depth. 
Figure  A-115  Lagrangian  Coordinates  When  Projectile  has  Penetrated  to  Depth  of 
4.4  m.  Radial  Dimensions  Doubled  for  Clarity. 

Figure  A-116  Normal  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  4.4  m Depth. 

Figure  A-117  Tangential  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  4.4  m Depth. 

Figure  A-118  Contours  of  Axial  Stress  in  Target.  Projectile  at  4.4  m Depth. 
Figure  A-119  Contours  of  Radial  Stress  in  Target.  Projectile  at  4.4  m Depth. 
Figure  A-120  Contours  of  Hoop  Stress  in  Target.  Projectile  at  4.4  m Depth. 
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Figure  A-121 

Figure  A-122 
Figure  A-123 
Figure  A -124 

Figure  A-125 

Figure  A-126 

Figure  A-127 
Figure  A-128 
Figure  A- 129 
Figure  A-130 

Figure  A-131 
Figure  A-132 
Figure  A-133 

Figure  A- 134 

Figure  A-135 

Figure  A-I36 
Figure  A- 137 
Figure  A-I38 


Contours  of  Radial-Vertical  Shear  Stress  in  Target.  Projectile 
at  4.4  m Depth. 

Contours  of  Vertical  Velocity  in  Target.  Projectile  at  4.4  m Depth. 
Contours  of  Radial  Velocity  in  Target.  Projectile  at  4.4  m Depth. 
Lagrangian  Coordinates  When  Projectile  has  Penetrated  to  Depth 
of  5.2  m.  Radial  Dimensions  Doubled  for  Clarity. 

Normal  Stress  (pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  5.2  in  Depth. 

Tangential  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  5-2  m Depth. 

Contours  of  Axial  Stress  in  Target.  Projectile  at  5.2  m Depth. 
Contours  of  Radial  Stress  in  Target.  Projectile  at  5.2  m Depth. 
Contours  of  Hoop  Stress  in  Target.  Projectile  at  5, .2  m Depth. 
Contours  of  Radial-Vertical  Shear  Stress  in  Target.  Projectile 
at  5*2  m Depth. 

Contours  of  Vertical  Velocity  in  Target.  Projectile  at  5.2  m Depth. 
Contours  of  Radial  Velocity  in  Target.  Projectile  at  5.2  m Depth. 
Lagrangian  Coordinates  When  Projectile  has  Penetrated  to  Depth 
of  6.8  m.  Radial  Dimens xons  Doubled  for  Clarity. 

Normal  Stress  (Pa)  Along  Surface  of  Projectile  When  Projectile 
is  at  6.8  m Depth. 

Tangential  Stress  (Pa)  Along  Surface  of  Projectile  when  Projectile 
is  at  6.8  m Depth. 

Contours  of  Axial  Stress  in  Target.  Projectile  at  6.8  m Depth. 
Contours  of  Radial  Stress  in  Target.  Projectile  at  6.8  m Depth. 
Contours  of  Hoop  Stress  in  Target.  Projectile  at  6.8  m Depth. 


Figure  A- 139  Contours  of  Radial-Vertical  Shear  Stress  in  Target.  Projectile 
at  6.8  m Depth. 

Figure  A-140  Contours  of  Vertical  Velocity  in  Target.  Projectile  at  6.8  m Depth. 
Figure  A-l4l  Contours  of  Radial  Velocity  in  Target.  Projectile  at  6.8  m Depth. 
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Figure  A-3 

Projectile  Depth  vs.  Time 
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Axial  Displacement  (m)  versus  Time  (s)  for  particle  at  Point  3, 
Z = 3-6  m,  R = 0.15  m. 
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Figure  A-19 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  4, 
Z ~ 0.4  m,  R = 0.15  m. 
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Figure  A-35 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  7,  Z - 5-4  m,  E = 0.45  m. 
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Figure  A-43 

Deviator  Stress,  J2,  History  (Pa)  for  Particle  ut  Point  8,  Z = 3.6  m,  R =,0.1(5  m 
(Octahedral  Shear  Stress  = v£/3(jP)) . 
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Figure  A-54 

Axial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  10,  Z = 0.4  m,  R = 0.45  m. 
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Figure  A-58 

)eviator  Stress,  J2,  History  (pa)  for  Particle  at  Point  10,  Z = 0.4  m,  R = 0.45  m 
.Octahedral  Shear  Stress  = V§/3(j2)  ) . 
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Figure  A-6l 


Radial  Displacement  (m)  versus  Time  (s)  for  Particle  at  Point  11,'  Z = 0,1  m,  R ■■  0.4>  m. 
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Figure  A-6i? 

Radial  Velocity  (m/s)  versus  Time  (s)  for  Particle  at  Point  12,-  Z = 3-6  m,  R *»  0.9  m. 
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POINT  12  TIME 

Figure  A-66 


Axial  Velocity  (m/e)  versus  Time  (s)  for  Particle  at  Point  12,  Z =3*6  m,  K “ 0.9 
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Figure  A-'fb 

Radial  Stress  (Pa)  History  for  'Particle  at  Point  13,  Z = 5.4  m,  R = 1,5  ffl, 
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Figure  A-75 

Pressure  (Pa)  versus  Time  (s)  for  Particle  at  Point  l4,  £ = 3.6  m,  R = 1.5  m. 
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Figure  A-85 

Contours  of  Radial -Vertical  Shear  Stress  in  Target.  Projectile  at  0.1  m Depth. 
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Figure  A-13i 

Contours  of  Vertical  Velocity  in  Target.  Projectile  at  5.2  m Depth. 
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Figure  A- 13?. 

Contours  of  Radial  Velocity  in  Target,  Projectile  at  5.2  m Depth 
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Figure  A-133 

Lagrangian  Coordinates  When  Projectile  has  Penetrated  to  Depth  of  6.8  m.  Radial 
Dimensions  Doubled  for  Clarity. 
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B-15.  Lagrangian  Coordinates  at  2.4  ms  (Radial  Dimensions  Doubled  to  Show 
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B-l6.  Contours  of  Radial  Velocity  at  2.4  ms  (Deformable  Penetrator  not 
Shown) . 
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B-80.  Vertical  Displacement  History  (m)  in  soil  at  Range  0,2'8  m,  Depth 
0.1  m (Deformable  Penetrator).  % 

B-8l.  Pressure  History  (Pa)  in  Soil  at  Range  0.28  m,  Depth  0.1  m (Deformable 
Penetrator) . 

B-82.  Deviator  Stress,  J2,  History  (Pa)  in  Soil  at  Range  0.28  m,  Depth 
0.1  m (Deformable  Penetrator). 

B-83.  Radial  Velocity  History  (m/s)  in  Soil  at  Range  0.28  ra,  Depth  0.2  m 
(Deformable  Penetrator). 
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B-84.  Vertical  Velocity  History  (m/s)  in  Soil  at  Range  0.28  m,  Depth 
0.2  m (Deformable  Penetrator) . 

B-85.  Radial  Displacement  History  (m)  in  Soil  at  Range  0.28  m,  Depth  0.2  m 
(Deformable  Penetrator) . 

B-86.  Vertical  Displacement  History  (m)  in  Soil  at  Range  0.28  m,  Depth  0.2  m 
(Deformable  Penetrator) . 

B-87.  Pressure  History  (Pa)  in  Soil  at  Range  0.28  m,  Depth  0.2  m (Deformable 
Penetrator) . 

B-88.  Deviator  Stress,  J2,  History  (pa)  in  Soil  at  Range  0.28  m,,  Depth 
0.2  m (Deformable  Penetrator). 

B-89.  Radial  Velocity  History  (m/s)  in  Soil  at  Range  0.28  m,  Depth  0.4  m 
(Deformable  Penetrator) . 

B-90.  Vertical  Velocity  History  (m/s)  in  Soil  at  Range  2.28  m.  Depth  0.4  m 
(Deformable  Penetrator) . 

B-91.  Radial  Displacement  History  (m)  in  Soil  at  Range  0.28  m,  Depth  0.4  m 
(Deformable  Penetrator) . 

B-92,  Vertical  Displacement  History  (m)  in  Soil  at  Range  0.28  m,  Depth  0.4  m 
(Deformable  Penetrator). 

B-93.  Pressure  History  (Pa)  in  Soil  at  Range  0.28  m,  Depth  0.4  m (Deformable 
Penetrator) . 

B-94.  Deviator  Stress,  J2,  History  (Pa)  in  Soil  at  Range  0.28  m,  Depth 
0.4  m (Deformable  Penetrator). 

B-95.  Radial  Velocity  History  (m/s)  in  Soil  at  Range  0.44  m,  Depth  0.1m 
(Deformable  Penetrator) . 

B-96.  Vertical  Velocity  History  (m/s)  in  Soil  at  Range  0.44  m,  Depth  0.1  m 
(Deformable  Penetrator) . 
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B-9'7.  Radial  Displacement  History  (m)  in  Soil  at  Range  0.44  m,  Depth  0.1  m 
(Deformable  Pcnetrator) . 

B-98.  Vertical  Displacement  History  (m)  in  Soil  at  Range  0.44  m,  Depth 
0.1  m (Deformable  Penetrator). 

B-99*  Pressure  History  (Pa)  in  Soil  at  Range  0.44  m,  Depth  0.1  m (Deformable 
Penetrator)  . 

B-100.  Deviator  Stress,  J2 , History  (Pa)  in  Soil  at  Range  0.44  m,  Depth 
0.1  m (Deformable  Penetrator). 

B-301.  Radial  Velocity  History  (m/s)  in  Soil  at  Range  0.44  in,  Depth  0.2  m 
(Deformable  Penetrator) . 

B-102.  Vertical  Velocity  History  (m/s)  in  Soil  at  Range  0.44  m,  Depth  0.2  m 
(Deformable  Penetrator). 

B-103.  Radial  Displacement  History  (m)  in  Soil  at  Range  0.44  m,  Depth  0.2  m 
(Deformable  Penetrator) . 

B-104.  Vertical  Displacement  History  (m)  in  Soil  at  Range  0.44  m,  Depth 
0.2  m (Deformable  Penetrator). 

B-105.  Pressure  History  (Pa)  in  Soil  at  Range  0.44  m,  Depth  0.2  m (Deformable 
Penetrator) . 


B-106.  Deviator  Stress,  J2,  History  (Pa)  in  Soil  at  Range  0.44  m,  Depth 
0,1  m (Deformable  Penetrator). 

B-107.  Radial  Velocity  History  (m/s)  in  Soil  at  Range  0.44  m,  Depth  0.4  m 
(Deformable  Penetrator). 

B-108.  Vertical  Velocity  History  (m/s)  in  Soil  at  Range  0.44  m,  Depth  0.4  m 
(Deformable  Penetrator) . 

B-109.  Radial  Displacement  History  (m)  in  Soil  at  Range  0.44  m,  Depth  0.4  m 
(Deformable  Penetrator) . 
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3-110.  Vertical  Displacement  History  (m)  in  Soil  at  Range  0.44  m,  Depth  0.4  m 
(Deformable  Penetrator) . 

b-111.  Pressure  History  (Pa)  in  Soil  at  Range  0.44  m,  Depth  0.4  m (Deformable 
Penetrator) . 

B-112 . Deviator  Stress,  J2,  History  (Pa)  in  Soil  at  Range  0.44  rr.,  Depth 

0.4  m (Deformable  Penetrator),  (octahedral  Shear  Stress  = v£73(J2)). 
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figure  B-ll 

Contours  of  radial  stress  at  1.1  ms  (deformable 


penetraior  not  shown). 
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B-12 

ile  penetrator  not  shown). 
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( Figure  B-22 

Deceleration  history  (g)  for  point  1 at  accelerometer  (filtered  above 

10  kH). 
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Figure  B-26 

Vertical  velocity  history  (m/s)  For  point  2,0.58  m behind  penetrator  tip  and  10  mm 
from  central  axis. 
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Figure  B-28 

DeceJ eration  history  (g)  for  point  5.28  ram  behind  pcnetrator  tip  and  6 ram  from 


centi  al  axi 


Figure  B-32 

history  (Pa)  for  point  3.0.79  m behind  penetrator  tip  and  69  ram 
Octahedral  shear  stress  = \j2/3  (J2)). 
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Figure  B-37 

Normal  stress  history  (Pa)  for  point  7,  0.l4  m behind  penetrator  tip  and  Ul  mm  from 
central  axis. 
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figure  B-^2 

Vertical  velocity  history  (m/s)  in  soil  at  range  75  min,  depth  0. 
penetrator) . 
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Radial  displacement  history  (m)  in  soil 
penetrator) . 
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Figure  B-W+ 

Vertical  displacement  history  (m)  in  soil  at  range  75  nun,  depth  0.1  m (de?ormable 
penetrator) . 
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•Figure  B-56 

Vertical  displacement  history  (m)  in  soil  at  range  75  nun,  depth  0.4  m (deformable 
penetrator) . 
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Figure  B-£5 

Radial  velocity  history  (m/s)  in  soil  at  ->”inge  O.lU  in,  depth  0.2  m (deformable 
penetrator) . 
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DEFORMABLE  BODY  FOR  WATCHING  HILL 


Figure  B-71 

Radial  velocity  history  (m/s)  in  soil  at  range  0.14  m,  depth  0.4  m (deformable 
penetrator ) . 
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j Figure  B-72 

j Vertical  velocity  history  (m/s)  in  soil  at  range  0.14  m,  depth  0.4  m (deformable 

penetrator) . 
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Figure  B 

Deviator  stress,  J2,  history  (Pa)  in  soil  at 
penetrator) , 
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Vertical  displacement  history  (m)  in 
penetrator) . 
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gure  B-86 

■oil  at  range  0.28  in,  depth  0.2  m (deformable 


Figure  B-89 

Radial  velocity  history  (m/s)  in  soil  at  range  0.28  m,  depth  0.4  m (deformable 
penetrator) . 
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Figure  B-96 

Vertical  velocity  history  (m/s)  in  soil  at  range  0.44  m,  depth  0,1  m (deformable 
penetrator) . 
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Figure  B-112 

Deviator  stress,  J2,  history  (Pa)  in  soil  at  rang; 
penetrator),  (octahedral  shear  stress  = ^ 2/3  (J2 
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